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Abstract 

Background: Breast mill< contains complex nutrients and facilitates the maturation of various biological systems in 
infants. Exosomes, membranous vesicles of endocytic origin found in different body fluids such as milk, can 
mediate intercellular communication. We hypothesized that microRNAs (miRNAs), a class of non-coding small RNAs 
of 18-25 nt which are known to be packaged in exosomes of human, bovine and porcine milk, may play important 
roles in the development of piglets. 

Results: In this study, exosomes of approximately 100 nm in diameter were isolated from porcine milk through 
serial centrifugation and ultracentrifugation procedures. Total RNA was extracted from exosomes, and 5S ribosomal 
RNA was found to be the major RNA component. Solexa sequencing showed a total of 491 miRNAs, including 176 
known miRNAs and 315 novel mature miRNAs (representing 366 pre-miRNAs), which were distributed among 30 
clusters and 35 families, and two predicted novel miRNAs were verified targeting 3'UTR of IGF-IR by luciferase assay. 
Interestingly, we observed that three miRNAs (ssc-let-7e, ssc-miR-27a, and ssc-miR-30a) could be generated from 
miRNA-offset RNAs (moRNAs). The top 10 miRNAs accounted for 74.5% (67,154 counts) of total counts, which were 
predicted to target 2,333 genes by RNAhybrid software. Gene Ontology and Kyoto Encyclopedia of Genes and 
Genomes (KEGG) pathway analyses using DAVID bioinformatics resources indicated that the identified miRNAs 
targeted genes enriched in transcription, immunity and metabolism processes, and 14 of the top 20 miRNAs 
possibly participate in regulation of the IgA immune network. 

Conclusions: Our findings suggest that porcine milk exosomes contain a large number of miRNAs, which potentially 
play an important role in information transfer from sow milk to piglets. The predicted miRNAs of porcine milk 
exosomes in this study provide a basis for future biochemical and biophysical function studies. 
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Background 

Milk, as the sole source of nutrition for infants, contains a 
potent mixture of diverse components such as milk fat 
globules (MFG) [1], immune competent cells and soluble 
proteins, for instance IgA, cytokines and antimicrobial pep- 
tides [2], which can provide protection against infections in 
newborns [3]. In addition, breast milk may have a role in 
tolerance induction [1] and may protect infants from devel- 
oping allergies [4]. 
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Exosomes are small (30-100 nm) membrane vesicles 
of endocytic origin that are released into the extracellu- 
lar environment upon fusion of multivesicular bodies 
(MVB) with the plasma membrane [5]. Many cells have 
the capacity to release exosomes, including reticulocytes 
[6], dendritic cells [7], B cells [8], T cells [9], mast 
cells [10], epithelial cells [11] and tumor cells [12]. In 
addition, exosomes have been found in physiological 
fluids, such as saliva [13,14], plasma [15], urine [16], am- 
niotic fluid [17], malignant ascites [18], bronchoalveolar 
lavage fluid [19] and synovial fluids [20]. Several studies 
have suggested that exosomes, which contain proteins, 
mRNA and microRNA (miRNA), stimulate and transfer 
surface receptors to target cells [21-23], as well as serve 
as novel vehicles for genetic exchange between cells 
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[24]. As with other biological fluids, microvesicle-like 
particles are also present in mouse milk [25] and human 
milk [26]. Recent published studies have isolated mRNAs 
and miRNAs fi:om bovine milk-derived microvesicles [27]. 
One study via deep sequencing technology identified 602 
unique miRNAs originating fi:om 452 miRNA precursors 
(pre-miRNAs) in human breast milk exosomes and found 
that, out of 87 well-characterized immune-related pre- 
miRNAs, 59 (67.82%) were enriched in breast milk exo- 
somes [28]. Recently, porcine milk was reported to contain 
180 pre-miRNAs, including 140 known and 40 novel por- 
cine pre-miRNAs, altogether encoding 237 mature miR- 
NAs [29]. 

MiRNAs are widespread among eukaryotes and represent 
key components of a conserved system of RNA-based gene 
regulation [30-33]. Many studies have demonstrated that 
miRNAs are key post-transcriptional regulators of gene 
expression and play important roles in a wide range of 
physiological and pathological processes [34], including 
development, differentiation, proliferation and immune re- 
sponses. It is believed that about 60% of mammalian genes 
are regulated by miRNAs [35-39]. 

Aside from being important farm livestock, pigs are 
also model animals for medical research. In the present 
study, we investigated miRNAs in milk exosomes of 
Landrace pigs in order to provide new information for 
investigations into the physiological functions of porcine 
milk. 

Methods 

Milk samples 

Porcine milk samples were collected between day 1 to 5 
after parturition from healthy lactating Landrace female 
pigs bred in the breeding farm of the Livestock Research 
Institute (Guangzhou, China). Milk samples were frozen 
immediately after milking and were kept at-80°C until 
use. 

Preparation of exosomes from milk 

Porcine milk samples were centrifiiged first at 2,000 x g for 
30 min at 4°C to remove MFGs as well as mammary gland- 
derived cells. Defatted samples were then subjected to 
centrifiigations at 4°C for 30 min at 12,000 x ^ to remove 
residual MFGs, casein and other debris. Subsequently, from 
the final supernatant (so-called whey or milk serum), the 
membrane fraction was prepared by ultracentrifiagation at 
110,000 x^ for 2 h in an SW41T rotor (Beckman Coulter 
Instruments, FuUerton, CA) [40]. 

Transmission electron microscopy (TEM) 

The final fraction obtained as described above was di- 
luted with 0.01 M PBS and ultracentrifuged again to 
recover microvesicles as pellets. Following fixation in 
2% glutaraldehyde, microvesicles were negatively stained 



with uranyl acetate and observed by TEM (JEOL 
JEM2000EX, Tokyo, Japan). 

RNA isolation and Solexa sequencing 

Total RNA was isolated from samples collected after 
ultracentrifugation using Trizol reagent (Invitrogen, 
Carlsbad, CA) according to the manufacturer s protocol. 
The quality of RNA was examined by 2% agarose gel 
electrophoresis and with a Biophotometer 6131 (Eppendorf, 
Germany), as well as further confirmed by using a 
Bioanalyzer (Agilent Technologies, Santa Clara, CA). 
Small RNAs (18-30 nt) were obtained from the total 
RNA, 5' and 3' adaptors were ligated to the small 
RNAs, then the adaptor-ligated RNAs were subse- 
quently transcribed into cDNA by RT-PCR, and the 
samples were amplified by PCR using primers comple- 
mentary to the two adaptors. The PCR products were 
purified and subjected to Solexa sequencing (Illumina, 
CA) at the Beijing Genomics Institute (BGI, Shenzhen, 
China). 

Sequence data analysis 

The raw reads obtained from Solexa sequencing were 
processed to obtain clean reads by summarizing data 
production, evaluating sequencing quality, calculating 
the length distribution of small RNA reads, removing 
low quality reads and adaptor sequences as described in 
previous paper [41]. All the clean reads were aligned 
against non-coding RNAs from the GenBank and Rfam 
(11.0) (ftp.sanger.ac.uk/pub/databases/Rfam) database to 
annotate and classify rRNA, tRNA, snRNA and other 
ncRNA sequences using tag2 annotation software (devel- 
oped by BGI). Then the selected sequences were mapped 
to the pig genome (sscrofa9, www.ensembl.org/Sus_scrofa/) 
using SOAPvl.ll software [42] to analyze their expression 
and distribution. Subsequently, the miRNA candidates were 
further analyzed by miRDeep 2 against all known miRNAs 
and porcine miRNA precursors (miRBase 20.0). All re- 
maining candidates who did not map to any miRNAs in 
miRBase 20.0 were considered as potential novel miRNAs. 
To further identify these potential novel miRNA candi- 
dates, software MIREAPvO.2 (http://sourceforge.net/pro- 
jects/mireap) [43] developed by BGI was used to predict 
novel miRNA by exploring the secondary structure, the 
Dicer cleavage site and the minimum free energy of the 
annotated small RNAs which could be mapped to genome. 
In briefly, the sequence length should be between 18-26 
nt, maximal free energy allowed for a miRNA precursor 
was -18 kcal/mol, maximal space between miRNA and 
miRNA* was 35 nt, and flank sequence length of 
miRNA precursor should be 10 nt. Finally, all remaining 
novel miRNA candidates were further subjected to 
MiPred (http://www.bioinf.seu.edu.cn/miRNA/) to filter 
out pseudo-pre-miRNAs. The minimum free energy 
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must be> -20 kcal/mol or P-value was >0.05 [44], and 
their secondary structures were also checked using the 
Mfold3.2 software [45]. All data for analysis in this study 
have been deposited in https://mynotebook.labarchives. 
com/share/allinchen/MTkuNXwxMzMxMS8xNS0yLlR 
yZWVOb2RlLzElNzEyODU2fDQ5LjU= with a DOI:10. 
6070/H4DN432G. 

PGR and qRT-PCR identification of known and novel miRNAs 

Total RNA (identical sample to that of the Solexa 
sequencing sample) was first digested with DNase I 
(Invitrogen), and 2 [ig of total RNA was reverse transcribed 
to poly (A) tail-added cDNA using the One Step Prime- 
Script miRNA cDNA Synthesis Kit (TaKaRa, Dalian) ac- 
cording to the manufacturers instructions. Briefly, a 60-nt 
adaptor containing a poly (A) structure was added to the 3' 
sequence of miRNAs, which were then reverse transcribed 
to an 80-bp cDNA sequence [46]. The cDNA was diluted 
5-fold with ddH20, and PGR was performed on a Bio-Rad 
system (BIO-RAD,USA )in a final 20 [iL volume reaction 
containing 2 \i\ PGR cDNA, 10 \iL of 2x PGR Mix (Toyobo, 
Osaka, Japan) and 1 mM of each primer. The real-time 
PGR thermal profile was as follows: 5 min at 95°G, 40 cycles 
of 30 s at 94°G, 30 s at the corresponding annealing 
temperature (Tm) and 72°G for 30 s, followed by 72°G 
at 10 min. PGR products were examined on an agar- 
ose gel to confirm that a single PGR product was gen- 
erated, and 5S ribosomal RNA was used as an internal 
control for the PGR. The miRNA forward primer was 
designed with Primer 5.0 (Table 1), and the reverse 
primer for miRNAs was the Uni-miR qPGR Primer of- 
fered by the kit One Step Prime Script miRNA cDNA 
Synthesis Kit (TaKaRa, Dalian). 

miRNAs target prediction and plasmid construction 

Two predicted novel miRNAs, named miR-PG-86 and 
miR-PG-263, were selected to predict their target genes 
in pig genome using the RNAhybrid software (http:// 
bibiserv.techfak.uni-bielefeld.de/rnahybrid/) with its own 
algorithm. The 3'-UTR sequences of porcine transcripts in 
whole genome were obtained from ensemble gene 66 
(sscorfa 9, www.ensembl.org/Sus_scrofa/). The 3'-UTR of 
IGF-IR contains the highly conserved binding sites for the 
two miRNAs, and the sequence (104 bp) is as follows: 
TCCTGGATCCCGATCCCGTGCAAACAGTACCGTGCG 
CACGCGGGCGGGCGGGGGGAGAGTTTTAACAATCT 
ATTCACAAGCCTCCTGTACCTCAGTGGATCTTG Fur- 
ther, the 3'-UTR sequence was inserted into pmirGLO 
Vector (Promega) with Xhol and Xbal double digestion to 
construct recombinant Dual-Luciferase reporter vector, 
named as pGLO-IGF-lR-3'UTR (Figure lA). Meanwhile, a 
plasmid containing mutant IGF-IR 3'-UTR, named as 
pGLO-IGF-lR-3'UTR-delete (Figure lA), was generated by 
deleting the core sequence of the two miRNA binding sites 



Table 1 PCR primers for miRNAs 



miRNAs name 


Primer sequence 


Renaturation 
temperature 


ssc-let-7e 


TGAGGTAGGAGGTOTATAG^ 


59.5°C 


ssc-miR-21 


GCTAGCTOTCAGACTGATGTO 


59.5°C 


ssc-miR-206 


TGGAATGTAAGGAAGTGTGTG 


59.5°C 


ssc- let-7i 


GCCGCTGAGGTAGTAGmGTGCT 


59.5°C 


ssc-miR-140 


GACAGTGG^ACCCTATGGTA 


59.5°C 


ssc-miR-92b-5p 


™tagggacgggacgcggtg 


59.5°C 


ssc-miR-22b-3p 


AAGCTGCCAGTOAAGAACTG 


59.5°C 


ssc-miR-28-5p 


GAAGGAGCTCACACTCTATOA 


59.5°C 


ssc-miR-205 


TCOTCATOCACCGGAGTCT 


59.5°C 


ssc-miR-451 


AAACCGTOCCATOCTGAG^ 


59.5°C 


ssc-miR-125b 


TCCCTGAGACCCTAAOTGTG 


59.5°C 


ssc-miR-9 


GCGGTCmGGTOTCTAGCTGT 


59.5°C 


ssc-let-7c 


TGAGGTAGTAGGTOTATGGT 


59.5°C 


P-m0281-5p(PS-16) 


TCTCCCAACCOTGTACCA 


58°C 


P-m0124-3p(PC-280) 


TGTOCGAGATOGGCTG 


58°C 


P-m0227-5p(PC-291) 


TOCTGAGTCGGACTGGG 


58°C 


P-m0355-5p(PC-82) 


CCCAGGATCAGAGGATGG 


58°C 


P-m0338-3p(PC-241) 


TCTGTGAACTAGAAACCTCTGG 


58°C 


P-m0105-3p(PC-72) 


CAmGATOAGTOGACACT 


58°C 


P-m0113-3p(PC-130) 


CTATGGATCTAGGAGGACGC 


58°C 


P-m0129-5p(PC-129) 


CTATGGATCTAAGAGGACACCC 


58°C 


P-m0058-5p(PC-276) 


TGTGTGTGATCGTOATGTGC 


58°C 


P-m0279-5p(PC-192) 


GTCOTGGTGAGTCGGATG 


58°C 


P-m0103-3p(PC-70) 


CATOCmGATCGTCTGG 


58°C 


P-m0265-3p(PC-139) 


CTGGAAGGAmGGGTAGG 


58°C 


P-m0210-5p(PC-277) 


TGTGTGTOTGTCGGATGAG 


58°C 


P-m0186-5p(PC-73) 


CA^AAGGATCGTGTGGG 


58°C 


P-m0070-3p(PC-63) 


CAGTAGGGATGAGAGGACACT 


58°C 



through DNA synthesis (Sangon Biotech (Shanghai) Go., 
Ltd.), the sequence is as follows: TCCTGGATCCCGATCC 
CGTGCAAACAGTACCGTGCGCACGCGGGCGGGCGG 
GGGGA GA GTTTTAACAATCTATTCACAA GCCTCCTG 
TACCC. 

Leuciferase reporter assay 

IPEG-J2 cells were maintained in DMEM/F12 (1:1) 
(GIBGO) and supplemented with 10% fetal bovine serum 
(FBS, GIBGO), 5 ng/ml EGF (peprotech, USA) and 5 ug/ 
ml insulin (Sigma, USA). Lipofectamine 2000 (Invitro- 
gen) was used for transfection. Gells (10,000) were plated 
in a 96-well plate. After 24 h cultivation, cells were 
transfected with a mixture including 500 ng pGLO-IGF- 
IR-3'UTR or pGLO-IGF-lR-3'UTR-delete construct and 
30pM of miR-PG-86 or miR-PG-263 mimics (GenePharma, 
Shanghai, Ghina). For control, 500 ng of pmirGLO- 
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TCAGTG GATCTT 



IGF-IR coding region 


3'UTR 


— IGF-IR mRNANM_214172 






Luciferase 


3'UTR 


— pGLO-IGF-lR-UTR 






Luciferase 


delete 


_ pGLO-IGF-lR-UTR -delete 



miR-PC-263 TGGGATCTGAGGACAGAAGGCGGA 
miR-PC-86 CCCATTGAGAGCGAAGGTGGCA 
IGF-IR mRNA(4133-4236) CTATTCACAAGCCTCCTGTACC TCAGTGGATCTT C 



IGF-lR-delete CTATTCACAAGCCTCCTGTACa 



C 3.50 

% 3.00 
re 

"5 2.50 

^ 2.00 

1 1.50 

re 1.00 
a) 

u 0.50 
0.00 

86mimics 263mimics NC 

pGLO-IGF-lR-3'UTR + + + 

iiiiR-pc-86 mimics + _ _ 




3.00 



> 2.50 




86mimics 263mimcs NC 

pGLO-IGF-lR-3'UTR-delete + + + 

miR-pc-86 mimics + — _ 



miR-pc-263 mimics 



miR-pc-263 mimics 



NC - - + NC - - + 

Figure 1 IVliR-PC-86 and IVIiR-PC-263 directly regulate IGF-IR expression via 3' UTR sites. (A) Schematic of IGF-IR mRNA and the luciferase 
reporter plasmids containing the miR-PC-86 and miR-PC-263 binding sites of IGF-IR mRNA. The 3' UTR sites were inserted downstream of the 
luciferase reporter, as indicated. TCAGTGG was the predicted target site of miR-PC-86, GGATC^ was the predicted target site of miR-PC-263. (B) 
miR-PC-86 and miR-PC-263 sequences and predicted binding site between miR-PC-86 and miR-PC-263 and IGF-IR mRNA. IGF-IR mRNA has one 
putative binding site for miR-PC-86/ miR-PC-263 on the 3' UTR. Twelve nucleotides TCAGTGGATC^ of IGF-IR 3' UTR (underlined) were delete in 
order to disrupt the binding with miR-PC-86 and miR-PC-263 seed regions. (C) IPEC-J2 cells were transfected with each of the constructed 
plasmids, together with miR-PC-86/ miR-PC-263and Renilla luciferase reporter plasmid (*P< 0.05, n = 6). 



scramble including a scrambled sequence of the miRNA 
target sequence was used. Cells were collected 48 h 
after transfection, and luciferase activity was measured 
using a Dual-GLO luciferase reporter assay system 
(Promega). Statistical differences between treatment 
and control groups were determined using Students 
t-test, at P < 0.05. 



Bioinformatics analysis 

Chromosomal localization and cluster analysis of miRNAs 

Pre-miRNAs of all miRNAs (known miRNAs and 
novel miRNAs) were mapped to the porcine genome 
(sscrofa9, www.ensembl.org/Sus_scrofa/) according to their 



positions on the chromosomes. Pre-miRNA positions less 
than 10 kb apart were considered to belong to the same 
miRNA cluster. 

Target prediction and Gene Ontology (GO) and Kyoto 
Encyclopedia of Genes and Genomes (KEGG) pathway 
analyses 

Porcine miRNA targets were obtained at the genome 
level. In brief, miRNA targets were predicted using the 
RNAhybrid software algorithm (http://bibiserv.techfak. 
uni-bielefeld.de/rnahybrid/) in 3'-UTR sequences of 
transcripts from the whole pig genome obtained from 
Ensembl Gene 66 database (sscrofa9, www.ensembl.org/ 
Sus_scrofa/). Strict criteria (perfect match of 2-8 nt in 
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Figure 2 Exosomes detected by TEM. The exosomes appeared as round or oval microvesicles (A, B), with a diameter of 50-150 nm and 
heavier density at the center than on the margin. 



the seed region; no more than -25 kcal/mol low free 
energy of miRNAs-mRNA binding) were appUed to the 
target prediction procedure. GO and KEGG pathway 
analyses were performed using DAVID bioinformatics 
resources (http://davidabcc.ncifcrf.gov/). 

Results 

Identification of exosomes 

Exosomes were obtained from porcine milk by ultra- 
centrifugation. After negative staining, approximately 
round-shaped porcine milk exosomes with diameters 
of -50-100 nm were observed by TEM, showing a 
greater density at the center than at the boundary 
(Figure 2A, B). 

Porcine milk exosomes contain RNA 

To further ascertain whether the porcine milk exosomes 
collected by ultracentrifugation contained RNA, we ex- 
tracted the samples using Trizol reagent and then exam- 
ined the recovered product by electrophoresis on a 2% 
agarose gel. To exclude the possibility of DNA contam- 
ination, total RNA was incubated at 37°C for 30 min 



with DNase I or RNase. The extracted RNA could not 
be digested by DNase I, while it could be degraded by 
RNase (Figure 3A). These results demonstrated that 
the porcine milk exosomes contained RNA. More inter- 
estingly, total RNA of porcine milk exosomes were 
enriched with 5S rRNA (Figure 3B), consistent with 
previous studies [27,28,47]. 

Solexa sequencing and analysis 
Solexa sequencing 

The sRNAs were enriched from porcine exosomes to con- 
struct a library for Solexa sequencing. We obtained 
9,033,167 raw reads and 6,013,724 high qualities reads after 
removal of low quality and contaminant reads. Among 
the high quality reads, there were 4,964,542 clean reads 
(82.55%), representing 1,691,655 unique sRNAs. The major- 
ity of the sRNAs in porcine milk exosomes were 18-25 nt 
in length (74.89%, Figure 4), with 2,458,894 reads (49.53%) 
representing 872,096 unique sRNAs (51.55%), including 
miRNAs and other sRNAs, such as rRNA, tRNA, snRNA, 
snoRNA, scRNA, small recognition particle RNA (srpRNA), 
repetitive sequence elements and unannotated sequences. 
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^ [FU]_ 
4- 




^ 5S 1 








2000bp — 






3- 






750bp — 
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lOObp — 






1- 






DNase 
RNase 






0- 






+ 

+ 




1 1 1 1 1 1 ' 
25 200 1000 4000 [nt] 


Figure 3 Milk-derived exosomes containing RNA. (A) Total RNA was extracted from porcine exosomes. M, 1, 2 and 3 represent the marker 
(DL 2000), RNA without any treatment, RNA treated with DNase 1 and RNA treated with RNase, respectively. (B) RNA sample analyzed by the 
Agilent Bioanalyzer 2100. 
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Figure 4 Length distribution of miRNAs reads from Solexa sequencing. A total of 4,964,542 clean reads were obtained, ranging from 10 to 
32 nt, most of which were 18-25 nt in length (accounting for 74.89%). 



which could be mapped to the pig genome. BLAST search- 
ing with the miRbase 20.0, identified a total of 230,216 reads, 
representing 1,555 unique known miRNAs. Due to RNA 
bias editing, 5' modifications and 3' modifications, many 
pre-miRNAs produce multiple mature miRNA isoforms, 
namely isomiRNAs, as described in many studies [48-50]. In 
the subsequent analysis, all isomiRNAs generated from the 
same precursor were considered one type of miRNA. Con- 
sequently, these 1,555 unique miRNAs corresponded to 176 
known mature miRNAs (205 pre-miRNAs, all the detail 
were listed in Additional file 1: Table SI). In addition, we 



identified 315 novel mature miRNAs (generated from 366 
pre-miRNAs, detail in Additional file 2: Table S2). Among 
the 315 novel miRNAs, 18 have not been deposited as por- 
cine miRNAs in miRbase 20.0, but had very high similarity 
with miRNA sequences of other species. These 18 miRNAs 
are labeled "PS" (Table 2), while 298 miRNAs that have not 
been deposited in miRbase 20.0 for any organism are labeled 
as "PC" and presented in Additional file 2: Table S2. 

There were 73 miRNAs with more than 100 counts 
and 264 miRNAs with less than 10 counts. The top 10, 
top 20, top 50 and top 100 miRNAs accounted for 



Table 2 Porcine novel miRNAs conserved in other species (miRBase release 20.0) 



Unique ID 


miRNAs name 


Count 


Sequence 


Size 


Conservation 


Match 


PS-1 


miR-290-5p 


14 


ACTCAAACTGTGGGGGCACm 


22 


mmu(#) 


Int sub(#) 


PS-2 


miR-378c 


33 


ACTGGAOTGGAGTCAGAAGT 


21 


hsa 


4nt delete 


PS-3 


miR-20b-3p 


4 


ACTGTAGTGTGGGCAOTCCAGT 


23 


hsa 


1 nt add, 1 nt sub 


PS-4 


miR-219-3p 


11 


AGAATOTGGCTGGACATCT 


20 


bta 


1 nt delete 


PS-6 


miR-138-5p 


42 


AGCTGGTGTOTGAATCAGGCCG 


23 


mmu 


perfect 


PS-7 


miR-31-5p 


5 


AGGCAAGATGCTGGCATAGCT 


21 


has 


perfect 


PS-9 


let-7f-l-3p 


2 


CTATACAATCTATOCOTCC 


21 


mo 


perfect 


PS-11 


miR-874-3p 


25 


CTGCCCTGGCCCGAGGGACCGA 


22 


mmu 


perfect 


PS-1 2 


miR-551a 


225 


GCGACCCACTOTGGmCC 


20 


hsa 


1 nt delete 


PS-1 3 


miR-138-3p 


1 


GCTAOTCACAACACCAGGGT 


21 


hsa 


Int sub, 1 delete 


PS-1 4 


miR-182-3p 


1 


GGTGGTOTAGAOTGCCAACT 


22 


mmu 


Int insert 


PS-1 5 


miR-5003-3p 


58 


TAmAATAGGTOTOGGA 


20 


hsa 


2nt sub, 2nt delete 


PS-1 6 


miR-150-5p 


164 


TCTCCCAACCOTGTACCAGT 


21 


mmu 


1 nt delete 


PS-1 7 


miR-2411-3p 


7 


TGAACTGTCATACTCCCACATC 


22 


bta 


3nt delete, 1 nt sub 


PS-1 8 


let-7f-5p 


20 


TGAGGTAGTAGATOTATAGTO 


23 


hsa 


1 nt insert 


PS-1 9 


miR-31-3p 


3 


TGCTATGCCAACATATOCCA 


21 


has 


1 nt delete 


PS-20 


miR-182 


46 


mGGCAATGGTAGAACTCACA 


22 


dre 


perfect 


PS-21 


miR-96-5p 


65 


mGGCACTAGCACAl II IIGCT 


23 


hsa 


perfect 



#: mmu, Mus musculus; hsa, Homo sapiens; bta. Bos taurus; dre, Danio rerio; rno, Rattus norvegicus; 

"sub", "delete", "add" represents nucleotide substitution, deletion, addition (at 5-end of miRNAs), respectively perfect stands for perfect match with reference miRNAs. 
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Figure 5 Distribution of miRNA reads and top 10 miRNAs. (A) The distribution of miRNA reads showed that the top 10, top 20, 

top 100 miRNAs accounted for 74.5%, 85.2% and 95.3% and 98.3% of total reads. (B) Cumulative proportions of top 10 miRNAs. MiR- 
ranked first, accounting for 29.6% of total reads. 



top 50 and 
-193-3p 



74.5%, 85.2%, 95.3% and 98.3% of the reads, respectively 
(Figure 5A). The top 10 miRNAs were ssc-miR-193a-3p, 
ssc-miR-423-5p, ssc-miR-320, ssc-miR-181a, ssc-miR-30a- 
3p, ssc-miR-378, ssc-miR-191, ssc-let-7a, ssc-let-7f and ssc- 
let-7c. With 67,154 counts (29.6%, average count: 460.6) 
(Figure 5B), ssc-miR-193a-3p ranked first among all 
miRNAs reads. 

Identification of miRNAs by PGR and sequencing 

To verify the deep sequencing results, we selected 
randomly 13 known miRNAs and 15 predicted novel 
miRNAs for PGR amplification (Figure 6A, B). Subse- 
quently, the 15 newly predicted miRNAs were cloned 
and sequenced. The results showed that 8 sequences 



5 6 7 8 9 10 11 12 13 14 




Figure 6 MiRNAs detected randomly in porcine milk. (A) Known 
miRNAs from miRBase (18.0), from M to 14, respectively: marker 
(DL 2000), ssc-let-7e, ssc-miR-21, ssc-miR-206, ssc-let-7i, ssc-miR-140, 
ssc-miR-92b-5p, ssc-miR-22-3p, ssc-miR-28-5p, ssc-miR-205, 
ssc-miR-451, ssc-miR-125b, ssc-miR-9, ssc-let-7c and 5 s (control). 
(B) Top 15 predicted novel miRNAs, from M to 16, respectively: 
marker (DL 2000), P-m0227-5p, P-m0338-3p, P-m0105-3p, 
P-m0058-5p, P-m0281-5p, P-m0265-3p, P-m0279-5p, P-m0103-3p, 
P-m0113-3p, P-m0129-5p, P-m0355-5p, P-m0210-5p, P-m0070-3p, 
P-m0124-3p, P-m0186-5p and 5 s (control). 



were fully matched, while 7 sequences had one or two 
mismatched nucleotides (Table 3). However their seed 
sequences remained unchanged. Simultaneously, the 
abundance of some novel miRNAs predicted by Solexa 
sequencing was confirmed by quantitative real-time 
PGR. The abundance of most miRNAs observed by 
qPGR of the sample pool and by sequencing were gen- 
erally consistent (Figure 7). 

Target verification of miR-PC-86 and miR-PC-263 against 
3'UTR of IGF-1R using luciferase report assay 

To investigate whether the predicted miR-PG-86 and 
miR-PG-263 (Figure 1) were functional novel miRNAs, 
target genes were predicted, and miR-PG-86/ miR-PG- 
263 were found to directly target IGF-IR 3'UTR se- 
quence. The full-length 3'UTR of IGF-IR mRNA was 
inserted downstream of the luciferase gene in the pmir- 
GLO Dual-Luciferase miRNA Target Expression Vector 
reporter plasmid, and the seed sequence was also delete 
to disrupt miR-PG-86/ miR-PG-263 binding (Figure IB). 
The wild-type (pGLO-IGF-lR-3'UTR) or delete (pGLO- 
IGF-lR-3'UTR-delete) plasmid was co-transfected with 
the miR-PG-86 and miR-PG-263 mimics into IPEG-J2 
cells. Forty-eight hours after transfection, the luciferase 
activity of the miR-PG-86 and miR-PG-263 group were 
significantly lower than that of the NG group (P < 0.05) 
respectively, and the reduction was rescued in the delete 
group (Figure IG). Thus, IGF-IR was initially confirmed 
as the target of miR-PG-86 and miR-PG-263. 

Genomic localization of pre-miRNAs 

To further establish the presence of miRNA precursors 
in the genome, all mature miRNAs (176 known and 315 
novels) were mapped to the S, scrofa genome (Figure 8). 
As a result, 176 known mature miRNAs were mapped to 
205 pre-miRNAs, and 315 novel miRNAs were mapped 
to 366 pre-miRNAs on the chromosomes. Our analysis 
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Table 3 mlRNAs matched to sequecing 

Predict new miRNA Matched sequence 

P-m028 1 -5 p T C T C C C AAC C C T T GT AC C ACT 

I I I I I I I I I I I I I I I I I I I I 
TCT CCCAACCCTTGT ACCAGT 

P-mOl 24-3p TGT TCCGAGATTGGG CTGCT 

I I I I I I I I I I I I I I I I I I 
TGTTCCGAGATTGGGCT GGC 

P-m0227-5p TTC CTGAGTCGGACT GGGCT 

I I I i I I I I I I i I I I I I I I I I 
TTC CTGAGTCGGACT GGGCT 

P-m0355-5p CCCAGGATCAGAGGATGGCTA 
I I I I I I I I I I I I I I I I I I I 
C C C AGGAT C AGAGGAT GGAC A 

P-m0338-3p TCTGTGAACTAGAAACCTCTGG 
I I I I I I I I I i I I i I i I I I I I I I 
TCTGTGAACTAGAAACCTCTGG 

P-mOl 05-3p CAT T T GAT TCAGTT GGAC ACT 

I I I I I I I I I I I I I I I I I I I I I 
CATTTGATTCAGTTGGACACT 

P-mOl 13-3p CTATGGATCTAGGAGGACGCA 
I I I I I I I I I I I I I I I I I I I I 
CTATGGATCTAGGAGGACGCC 

P-m0129-5p CTATGGATCTAAGAGGACACCC 
I I I I I I I I I I I I I I I I I I I I I I 
CTATGGATCTAAGAGGACACCC 

P-m0058-5p TGT GTGTGATCGTTAATGTGC 

I i I I I I I I I I I I I I I I I i I I I 
TGT GTGTGATCGTTAATGTGC 



Table 3 mIRNAs matched to sequecing (Continued) 

P-m0279-5p GTC CTTGGTGAGTCGGATGTC 

I I I I I I I I I I I I I I I I I I I I I 
GTCCTTGGTGAGTCGGATGTC 

P-mOl 03-3p CATTGCTTTGATCGTCTGGC 
I I I I I I I I I I I I I I I I I I I I 
CAT TGCTTTGATCGT CTGGC 

P-m0265-3p CTGGAAGGATTTGGGTAGGGA 
I I I I I I I I I I I I I I I I I I I 
CTGGAAGGATTTGGGTAGGTC 

P-m0210-5p TGTGTGTTCTGTCGGATGAGA 
I I I I I I I I I I I I I I I I I I I I 
TGT GTGTTCTGTCGGATGAGC 

P-mOl 86-5p CAT TTTAAGGATCGT GTGGGT 

I I I I I I I I I I I I I I I I I I I I 
CAT TTTAAGGATCGT GTGGGG 

P-m0070-3p CAGTAGGGATGAGAGGACACT 
I I I I I I I I I I I I I I I I I I I I I 
C AG T AGGGAT GAGAG GAC AC T 



revealed that the genomic density distribution of porcine 
milk pre-miRNAs (number of pre-miRNAs per Mb of 
each chromosome) was heterogeneous (Figure 8), ran- 
ging from 0.45 to 0.11 pre-miRNAs for 1 M of genomic 
sequence. Chromosomes with the highest and lowest 
densities of pre-miRNAs were chromosome 12 (29 pre- 
miRNAs per 64 Mbp) and chromosome 13 (25 pre- 
miRNAs per 219 Mbp), respectively. Interestingly, the 
medium-length X chromosome (144 Mbp, ranking 10th 
in length among the 19 chromosomes in pigs) was an 
exception by encoding an intermediate number (25 out 
of 366, 6.8%) of pre-miRNAs, corresponding to 0.17 pre- 
miRNAs for 1 M of genome sequence, but yet contained 
the most clusters. 

In addition, we observed many mature miRNAs having 
multiple miRNA precursors located in the same or dif- 
ferent chromosomes. Of the novel predicted miRNAs, 
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Figure 7 Expression of 15 predicted novel miRNAs in the sample pool detected by qRT-PCR. Trends in relative expression by qRT-PCR and 
counts from Solexa sequencing of miRNAs, except for PC-192, were consistent. 



40 pre-miRNAs had two copies in the genome, 7 pre- 
miRNAs had 3 copies, 2 pre-miRNAs had 4 copies, 1 
pre-miRNA had 8 copies and 249 pre-miRNAs were 
unique. With regard to known miRNAs, 4 pre-miRNAs 
had 3 copies, 22 pre-miRNAs had two copies and 149 
pre-miRNAs had only one copy. 

Mature miRNAs 

It is well accepted that only one of two strands gener- 
ated from a precursor is preferentially incorporated into 
RNA-induced silencing complexes (RISC), whereas the 



complementary strand (miR*) may be degraded. Closer 
examination of mature miRNAs generated from pre- 
miRNAs showed that precursor miRNAs could be di- 
vided into three groups (Table 4): pre-miRNAs only with 
the left- arm sequence (5p), pre-miRNAs only with the 
right-arm sequence or both. Most pre-miRNAs seemed 
to be single-arm miRNAs (5p or 3p), while 50 pre- 
miRNAs possessed both 5p and 3p sequences (40 coupled 
mature miRNAs, all details were listed in Additional file 3: 
Table S3). Further analysis of the 40 coupled mature 
miRNAs indicated that most of the pre-miRNAs had no 
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Figure 8 Distribution of 30 miRNA clusters. Tine number of base points near different bars indicates tine number of clusters in the chromosome. 
The relative vertical dimension of the point on the bar represents the location of cluster. The label "numberl /number2" above every bar indicates the 
value of "pre-miRNAs/mature miRNA". 
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Table 4 Pre-miRNAs and their corresponding mature 



miRNAs type 





Pre-miRNAs 


miRNA-5p 


miRNA-3p 


Both 


Mature 
miRNAs 


known 


205 


67 


57 


33(26)# 


176 


novel 


366 


150 


139 


17(13) 


315 


total 


571 


217 


196 


50(39) 


491 



#Number in bracl<et represents mature miRNAs couple number. 



significant difference in abundance between the 5p-arm 
and 3p-arm sequences (Figure 9). Some miRNAs had dif- 
ferent expression levels between the two strands. For ex- 
ample, ssc-miR-193a-3p had 67,154 counts, ranking first 
among all miRNAs, while ssc-miR-193a-5p had 2,538 
counts. Conversely, miR-423-5p had 22,588 counts, while 
its complementary strand, miR-423-3p, had only 654 reads, 
and which was shared by miR-22, miR-30a, miR-339, 
miR-17, miR-24, miR-331, miR-27b and let-7d (detail in 
Additional file 3: Table S3). 

As described in other studies [51], some small RNAs 
were generated from the loop or the region between the 
loop and stem (ssc-let-7e, ssc-miR-27a and ssc-miR-30a) 
(Figure lOA-C). In addition, we detected another interest- 
ing type of small RNAs known as miRNA-offset RNAs, or 
"moRNAs", which are derived from the ends of pre- 
miRNAs, predominantly from the 5' end, independent of 
the mature miRNA. A good example of moRNAs and 
small RNAs generated from the loop was pre-miR-30a 
(Figure IOC). At the 5' end of pre-miR-30, a 18 nt RNA 
sequence was found to be generated from the loop, down- 
stream of ssc-miR-30a-5p (Figure IOC). The findings sug- 
gest that these RNAs may be only byproducts of Drosha 
and Dicer processing, or these small RNAs alternatively 



may take part in other important regulatory functions dif- 
ferent from those of miRNAs. 

MiRNA clusters 

According to criteria for classification in miRbase, pre- 
miRNAs located on a chromosome with an interval of 
less than 10 kb are defined as belonging to an miRNA 
cluster. In our analysis, 30 (including 11 novel and 19 
known) (Figure 8) clusters were detected (Table 5). 
Among all miRNAs clusters, there were several pre- 
miRNAs with intervening sequences of less than 1 kb, 
including 10 known clusters (miR-99b/let-7e/125a, 
miR-24-2/27b/23b, miR-99a/let-7c, miR-29b/29a, miR-221/ 
222, miR-98/let-7f, miR-181c/d, miR-363/92a/19b-2/106a, 
miR-363/92a/19b-2, miR-181b-l/181a-l and miR-17/18a/ 
19b-l/92a-l) and 4 novel miRNAs clusters (cluster 3, 9, 12, 
22). We identified a typical polycistronic miRNA cluster, 
miR-363/92a/19b-2/20b, on chromosome X. Interestingly, 
the homologous cluster, miR-363/92a/19b-2/20b/106a on 
chromosome X, was located 33.5 kb downstream of miR- 
363/92a/19b-2/20b (Figure llA, C), and a paraologous clus- 
ter miR-17/18a/19b-l/92a-l was found on chromosome 11 
(Figure IIB, C). The organization of miRNA precursors in 
the genome may account for variable levels of expression 
and regulation of mature miRNAs. 

MiRNA families 

It is widely believed that the members of a given miRNA 
family regulate very similar sets of target genes. Apart 
from miRNA clusters, miRNA families were also recog- 
nized in the miRNAs of exosomes. Based on seed se- 
quences, 35 miRNA families were identified (26 known 
and 9 novel miRNA families) to contain at least two 
members and the identification of novel miRNAs added 
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Figure 9 5p and 3p arm expression of 46 pre-miRNAs. (A) 30 known pre-miRNAs. (B) 1 1 novel pre-miRNAs (representing 15 coupled 3p and 
5p arm sequences). 
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B 
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Figure 10 Three distinctive pre-miRNAs identified in porcine milic exosomes. (A) The ssc-miR-27a precursor produced a 3p-arm miRNA 
sequence and a loop-derived small RNA. (B) The ssc-let-7e precursor produced a 5p-arm sequence, a Sp-arm mlRNA sequence and a loop-derived 
small RNA. (C) The ssc-miR-30a precursor produced a 5p-arm sequence, a 3p-arm miRNA sequence and a loop-derived small RNA. In addition, 
ssc-moRNA-3, belonging to new type of miRNA termed moRNA, was found at the 5' end of pre-miR-30. 



new members to 5 known families (Table 6). In our 
study, 8 miRNA families (let-7, mir-1, mir-17, mir-181, 
mir-148, mir-30, mir-92 and mir-99) were found with at 
least 3 members among all exosome miRNAs. The let-7 
family had 9 members, miR-181 family had 4 members 
(miR-181a/b/c/d) and miR-30 family had 5 members 
(miR-30a/b/c/d/e). Most importantly, these miRNAs 
were highly expressed, and let-7 family members (let-7a/ 
f/c), miR-181a and miR-30a-3p were enriched among 
the top 10 miRNAs. However, members in the same 
family were highly differentially expressed. In the 
miR-181 family, miR-181a and miR-181b were dominant 
types with 13,345 reads and 3,333 reads, respectively. 
Similarly, miR-30a was the most abundant in the miR-30 
family. The differential expression of members in the 
same family may be partly due to regulation of their pre- 
cursors [52]. On the other hand, combined with the 
cluster analysis, we also observed that some miRNAs 
shared not only the same cluster but also the same fam- 
ilies. These miRNAs included 181a/b, let-7f/miR-98, 
181c/d, let-7a/f-5p/d-5p, 30b/d, 30c/e and miR-221/222. 
More interestingly, family members of the same cluster 
seemed to share expression patterns (Figure 12A-E). 
As mentioned above, miR-17-5p, miR-363, miR-106a, 



miR-18a, miR-19b, miR-92a, miR-20b and miR-92b formed 
a complex cluster and family network, and they also 
showed different expression patterns. MiR-92a, miR-19b 
and miR-363 were found to be highly expressed, while 
miR-17-5p, miR-18a, miR-20b and miR-106a were 
lowly expressed. The difference in abundance of the 
homologous or paralogous clusters may be attributed 
to the copy number of miRNA precursor itself or to 
the post-transcriptional regulation of the process 
of generating a mature miRNA from the precursor 
miRNA. 

In addition, many miRNA families showed low expres- 
sion (count number <100) in milk exosomes, such as the 
miR-1, miR-130, miR-17, miR-10, miR-29, miR-374, 
mir-9, miR-1 5 and miR-491 families (Figure 12F), which 
are routinely expressed in specific tissues [53-56]. Inter- 
estingly, all 9 novel miRNAs families showed extremely 
low expression levels (Figure 12G, count number <50), 
indicating that these miRNAs may only be expressed in 
certain physiology processes and may be the reason for 
why these miRNAs have not been detected until now. 
MiR-148a was reported to be an important biomarker 
for milk exosome miRNAs [28,57]. In this study, three 
members of this family (miR-148a, miR-148b and 
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Table 5 miRNAs cluster 



Table 5 miRNAs cluster (Continued) 



Cluster Cluster name 
No. 



members 



chromosome 



uster 1 new 

uster 2 mir-181a-2/181b-2 

uster 3 new 

uster 4 miR-24-l/27a/23a 



uster 5 

uster 6 

uster 7 

uster 8 

uster 9 
uster 10 

uster 1 1 

Cluster 12 
cluster 13 
cluster 14 

uster 15 

uster 16 



miR-181c/d 

miR-143/145 

let-7a/let-7f-2/let- 
7d 

miR-30d/30b 



mir-99b/let-7e/ 
125a 

mir-30c-l/ mir- 
30e 

new 

let-7a-l/miR-100 
new 

miR-181b-l/181a- 
1 

miR-24-2/27b/23b 



uster 17 miR-17-92a-l 

uster 18 new 

uster 19 miR-425/191 

uster 20 miR-99a/let-7c 

uster 21 new 

uster 22 new 

uster 23 new 

uster 24 new 

uster 25 miR-29b-l/29a 

uster 26 new 

uster 27 miR-221/222 



p-m0165, p-m0166 1 

ssc-mir-181a-2, 1 
ssc-mir-181 b-2 

ssc-miR-199b, P-m0179 1 

ssc-mir-24-1, 2 

ssc-miR-27a, 

ssc-mir-23a 

ssc-miR-181c, 2 
ssc-mir-181d 

ssc-mir-143, 2 
ssc-miR-145 

ssc-let-7a-2, P-m0204, 3 
ssc-let-7f-2, ssc-let-7d 

ssc-mir-30d, ssc-mir- 4 
30b 

P-m0263, P-m0264 6 

ssc-mir-99b, ssc-let-7e, 6 
ssc-mir-125a 

ssc-mir-30c-l, ssc-mir- 6 
30e 

P-m0302,P-m0318 7 

ssc-let-7a-l, ssc-mir-100 9 

P-m0346, P-r7i0347, 9 
P-m0348,P-m0349 

ssc-mir-181b-l, 10 
ssc-mir-181a-l 

ssc-mir-24-2, 10 

ssc-mir-27b, 

ssc-mir-23b 

ssc-mir-17, 11 
ssc-mir-18a, 
ssc-nnir-19b-l, 
ssc-nnir-92a-l 

P-m0024, P-m0025, 12 
P-m0026 

ssc-mir-425, 13 
ssc-mir-191 

ssc-mir-99a, 13 
ssc-let-7c 

P-m0099, 15 
P-m0083 

P-m0124, 17 
P-m0113 

P-m0125, P-m0126, 17 
P-m0127, P-r7i0128, 
P-m0114 

P-m0130, P-m0131 17 

ssc-mir-29b-l, ssc-mir- 18 
29a 

P-m0136, P-m0137, 18 
ssc-mir-183 

ssc-mir-221, ssc-mir-222 X 



cluster 28 


miR-98/let-7f-l 


ssc-mir-98, ssc-let-7f-l 


X 


cluster 29 


miR-363~20b 


ssc-mir-363-1, 
ssc-mir-92a-2, 

bbL 1 1 III 1 y\J Z, 

P-m0370, 

sc-mir-20b- 

l/ssc-mir-20b-2 


X 


Cluster 30 


miR-363~20b 


ssc-mir-363-1, 

ssc-mir-92a-2, 

ssc-mir-19b-2, 

P-m0371, 

ssc-mir-20b-l/ssc- 

mir-20b-2, 

ssc-miR-106a 


X 



miR-152, Figure 12H) showed modest expression levels, 
suggesting that miR-148 may be a stably expressed 
miRNA in exosomes of most mammals including pigs. 

GO and KEGG pathway analyses 

To better understand the role of miRNAs in pig milk 
exosomes, potential targets of miRNAs at the pig gen- 
ome level were explored by using our previous method 
[58]. We selected the top 10 miRNAs for target predic- 
tion using RNAhybrid, which found 2,333 potential tran- 
scripts. GO analysis showed that those targets were 
enriched in many processes (Table 7), including antigen 
processing and presentation, MHC class II protein com- 
plex and transcription. Pathways analysis revealed these 
targets were enriched in 7 KEGG pathways (Table 8), 
which are essential for piglet metabolism (glyceropho- 
spholipid metabolism and citrate cycle), immunity 
(cell adhesion molecules, asthma and intestinal immune 
network for IgA production, antigen processing and 
presentation) and development (Notch signaling path- 
way). The results suggest that these milk miRNAs likely 
take part in regulation of important protective functions 
in piglets, including the intestinal immune network for 
IgA production and antigen processing and presentation. 

Targets of miRNAs in IgA immune network 

As shown by GO and KEGG analysis, many of the iden- 
tified miRNAs were predicted to participate in immun- 
ity, similar to findings of other studies [28]. To further 
understand those miRNAs, the top 20 predicted 
miRNAs were selected for in-depth analysis of the IgA 
immune network. The results showed that 14 of the top 
20 miRNAs are likely involved in the IgA network and 
may target about 20 immune-related genes {APRIL, 
CCL2S, CD40, CD80, CD86, ICOS, ILIO, IL-S, ITGB7, 
LOC100621SS9, LTBR, MADCAMl, SLA, SLA-DMA, 
SLA-DMB, SLA'DOB, SLA-DQAl, SLA-DRA, SLA-DRBl 
and TGFB3) (Figure 13). These results suggest that por- 
cine milk miRNAs take part in regulation of the IgA 
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Figure 11 MiR-363~20b~106a homologous or paralogous cluster and expression level. (A) The miR-363/92a/19b/20b cluster and its 
homologous cluster miR-363/92a/19b/20b/106a on chromosome X were separated by a 33.5 kb DNA fragment. (B) The miR-17/18a/19/92a cluster 
was located on chromosome 1 1. In the genome, miR-92a/19b showed three copies; miR-363 and miR-20b had two copies; while miR-17, miR-18a 
and miR-106a had one copy. (C) Expression of mature miRNAs produced from the miR-363~20b~106a cluster and miR-17~92 cluster. 



network and immunity of piglets. In addition, several 
miRNAs shared the same target gene. Interestingly, all 
of the let-7 family members (let-7a, let-7c and let-7f) 
could target CCL2S, Those miRNAs are proposed to 
play a key role in IgA production in the piglet digestive 
tract and deserve further exploration, as mucosal im- 
munity is critically important for the protection of 
newborn piglets. 

Discussion 

In the present study, a comprehensive miRNA expres- 
sion profile of porcine breast milk exosomes was ex- 
plored via a deep sequencing approach. We found in 
total 176 known miRNAs (miRBase 20.0) and 366 pre- 
miRNAs producing 315 mature miRNAs. Luciferase 
reporter assay was used to explore the targets of two 
predicted novel miRNAs in this study. Results indicated 
both of them down-regulated the luciferase expression 
by targeting 3'UTR of IGF-IR. All these pre-miRNAs 
were distributed in 30 clusters (11 novel and 19 known 
clusters), and the mature miRNAs could be assigned to 
35 families (26 known and 9 unknown families). GO and 
KEGG pathway analyses show that those miRNAs may 
participate in many different immune-related processes. 
An analysis of the top 20 miRNAs showed that 14 of 
them may be involved in many regulatory aspects of the 
IgA immune network. 



A recent study of exosome miRNAs in Yorkshire sow 
milk discovered 180 pre-miRNAs, including 140 known 
porcine pre-miRNAs and 40 novel pre-miRNAs, which 
encode 237 mature miRNAs (234 unique miRNAs) [29] 
In the current study, we discovered 205 known porcine 
pre-miRNAs (176 mature miRNAs) and 366 novel pre- 
miRNAs (315 mature miRNAs), approximately 254 more 
mature miRNAs than were revealed in the former report 
[29]. Therefore, our results substantially supplement the 
known pig miRNAs, particularly milk exosome miRNAs. 
Interestingly, most of the novel miRNAs were low in 
abundance (312 miRNAs with less than 100 reads and 
only 3 miRNAs with >100 reads), which is possibly the 
reason for why these miRNAs were not detected in a 
previous study by Gu et al, [29]. Further comparison re- 
vealed that miR-191 and let-7a, which potentially play a 
vital role in immunity, were found in both that study 
and in the top 10 miRNAs of our study. Other miRNAs 
identified previously (miR-30a-5p, miR-25-3p, miR-182- 
5p, miR-200c-3p and miR-375-3p) were not detected in 
our study. Furthermore, miR-148a, a potential biomarker 
for quality control in bovine milk [57] and human milk 
[28], which was found to be highly expressed throughout 
the lactation period of Yorkshire sows [29], was only de- 
tected at a moderate level in Landrace pigs in our study. 
MiR-148a has been reported to be a tumor metastasis 
suppressor in gastric cancer [59], and ectopic expression 
of miR-148a was shown to induce apoptosis and silence 
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Table 6 mlRNAs family 





Family 


Seed 


Number 


Members 


Family 1 


let-7 


GAGGTAG 


9 


ssc-let-7i, ssc-miR-98, 
ssc-let-7a,ssc-let-7f, 
ssc-let-7c, ssc-let-7 g, 
ssc-let-7e, ssc-let7d-5p 


Family 2 


mir-1 


GGAATGT 


3 


ssc-miR-1, ssc-miR-206, 
PC-1 17 


Family 3 


mir-10 


ACCCTGT 


2 


ssc-miR-lOb, 
ssc-miR-lOa 


Family 4 


mir-1 03 


GCAGCAT 


2 


ssc-miR-107, 
ssc-miR-103 


Family 5 


mir-1 25 


CCCTGAG 


2 


ssc-miR-125a, 
ssc-miR-125b 


Family 6 


mir-1 30 


AGTGCAA 


2 


ssc-miR-130b, 
ssc-miR-130a 


Family 7 


mir-1 48 


CAGTGCA 


3 


ssc-miR-148b, 

ssc-miR-152, 

ssc-miR-148a 


Family 8 


mir-1 7 


AAAGTGC 


3 


ssc-miR-106a, 

ssc-miR-17-5p, 

ssc-mir-20b-l/ssc-mir-20b-2 


Family 9 


mir-1 81 


ACATOA 


4 


ssc-miR-181a, 
ssc-miR-181b, 

rrr- k-v-, ID 1 O 1 ^ 

ssc-miK- 1 O 1 C, 

ssc-miR-181d-5p 


Family 10 


mir-1 86 


AAAGAAT 


2 


ssc-miR-186, PC-36 


Family 1 1 


mir-1 93 


ACTGGCC 


2 


ssc-miR-193a-3p, PC-3 


Family 12 


mir-221 


GCTACAT 


2 


ssc-miR-222, ssc-miR-221 


Family 13 


mir-23 


TCACA^ 


2 


ssc-miR-23b, ssc-miR-23a 


Family 14 


mir-27 


TCACAGT 


2 


ssc-miR-27a,ssc-miR-27b 


Family 15 


mir-29 


AGCACCA 


2 


ssc-miR-29b, ssc-miR-29c 


Family 16 


mir-30 


GTAAACA 


5 


ssc-miR-30e-5p, 

ssc-miR-30c, 

ssc-miR-30d, 

ssc-miR-30b-5p, 

ssc-miR-30a-5p 


Family 17 


mir-30(#) 


mCAGT 


2 


ssc-miR-30a-3p, 
ssc-miR-30e-3p 


Family 18 


mir-339 


CCCTGTC 


2 


ssc-miR-339-5p, 
ssc-miR-4334-3p 


Family 19 


mir-34 


GGCAGTG 


2 


ssc-miR-34c, 
ssc-miR-34a, 


Family 20 


mir-363 


A^GCAC 


3 


ssc-miR-363, 

ssc-miR-92b-3p, 

ssc-miR-92a 


Family 21 


mir-374 


TATAATA 


2 


ssc-miR-374b-5p, 
<;<;r-miR--^74?i 


Family 22 


mir-378 


CTGGACT 


2 


ssc-miR-378, PS-2 


Family 23 


mir-491 


GTGGGGA 


2 


ssc-miR-491, PC-1 22 


Family 24 


mir-497 


AGCAGCA 


4 


ssc-miR-497, ssc-miR-15b, 
ssc-miR-16, ssc-miR-15a 


Family 25 


mir-9 


CmGGT 


2 


ssc-miR-9-2, ssc-miR-9-1, 


Family 26 


mir-99 


ACCCGTA 


3 


ssc-miR-99a, 
ssc-miR-99b, 
ssc-miR-100 



Table 6 mIRNAs family (Continued) 



Family 27 


new-1 


CATGA^ 


2 


P-m0040-3p(PC-9), 
P-m0240-5p(PC-232), 


Family 28 


new-2 


AGAGGGA 


2 


P-m0048-3p(PC-152), 
P-m0064-3p(PC-217) 


Family 29 


new-3 


AmCAT 


2 


P-m0084-5p(PC-228), 
P-m0105-3p(PC-72), 


Family 30 


new-4 


GCTAGGA 


2 


P-m0110-5p(PC-112), 
P-m0139-5p(PC-lll) 


Family 31 


new-5 


TATGGAT 


2 


P-m0113-3p(PC-130), 
P-mOl 29-5 p(PC- 129) 


Family 32 


new-6 


CCTGGAT 


2 


P-mOl 42-5 p(PC- 164), 
P-m0179-3p(PC-165) 


Family 33 


new-7 


CATAm 


2 


P-m0161-3p(PC-231), 
P-m0244-3p(PC-162) 


Family 34 


new-8 


GmCGA 


2 


P-m0183-3p(PC-286), 
P-m0325-5p(PC-128) 


Family 35 


new-9 


cmcGG 


2 


P-m0229-5p(PC-15), 
P-m0342-3p(PC-105) 



The underline indicated this family contains novel miRNAs. #: due to miRNAs 
classification by seed sequence, 3p and 5p of miR-30 represent different 
miRNAs families. PC is unique ID for porcine miRNAs candidate. 



Bcl-2 in colorectal cancer cells [60]. By bioinformatics 
analysis, miR-148a was determined to be possibly related 
to immunity and gastrointestinal health, but the under- 
lying regulatory mechanism remains unclear. 

MiR-92a belongs to the miR-17 -92 cluster with seven 
miRNAs (miR-17-5p, miR-17-3p, miR-18a, miR-19a, 
miR-19b, miR-20a and miR-92a) and was first described as 
an oncogenic miRNA cluster involved in B-cell lymphoma 
[61]. Recent studies indicated that the miRNA-17-92 
(miR- 17-92) cluster directly targets the TGFB pathway in 
cancer cell lines in the mouse embryo stage [62]. In 
addition, the miR- 17-92 cluster also participates in normal 
development of the heart, lungs and immune system [63]. 
MiR- 19 can promote leukemogenesis in Notchl -induced 
T-cell acute lymphoblastic leukemia (T-ALL) in vivo [64]. 
Overexpression of the mir-17-mir-18a-mir-19b-l cluster 
was shown to accelerate Myc-induced tumor development 
in a mouse B-cell lymphoma model [61]. The combined re- 
sults above imply that members of the cluster miR-363/ 
92a/ 19b-2/20b/ 106a may be related to cell proliferation and 
development. In porcine milk, miR-363/92a/19b-2/20b 
(miR-363/92a/19b-2/20b/106a) and miR-17/18a/19b-l/92- 
1 were also detected. The miR- 181 (181a/b/c/d) family is 
related to the development of different cells. It was reported 
that miR-181c/d can inhibit cell cycle and proliferation and 
that miR-181c regulates TNF-a [65]. The miR-30(b/c/d/e) 
family regulates kidney development by targeting the tran- 
scription factor Xliml/Lhxl in Xenopus [66]. The well- 
known let-7 (a/b/d/f) family is involved in oncogene expres- 
sion [67], and let-7/miR-98 family members are expressed 
late in mammalian embryonic development [68]. Thus, 
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Figure 12 Expression of miRNA families. (A) Expression of the I et-7f family; let-7f-l, let-7a and let-7d formed a cluster, while miR-98 and let-7e 
formed a cluster. (B) Expression of miR-30 family; miR-30d clustered with miR-30b, and miR-30c clustered with miR-30e. (C) Expression of miR-181 
family; miR-181 a clustered with miR-181b, and miR-181c clustered with miR-181d. (D) MiR-222 and miR-221; miR-222 and miR-221 belonged to 
the same family and the same cluster. (E) MiR-27 family and miR-23 family; miR-27b and miR-23b formed a cluster in the genome. (F) Expression 
and clustering of remaining miRNA families; members 1, 2, 3 and 4 represent members of the indicated family. (G) Novel miRNA family with low 
expression. (H) Expression of the miR-148 family containing milk "marker" miRNAs. 
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Table 7 Gene ontology analysis of potential targets of 



top10 mlRNAs 



Category 


Term 


Count 


P-Value 


Biological process 


transcription 


15 


6.70E-03 


Biological process 


regulation of transcription 


21 


1 .OOE-02 


Biological process 


antigen processing and 
presentation 


7 


1 .50E-02 


Biological process 


regulation of RNA 
metabolic process 


17 


1 .70E-02 


Biological process 


regulation of transcription, 
DNA-dependent 


17 


1 .70E-02 


Biological process 


antigen processing and 
presentation of peptide 
or polysaccharide antigen 
via MHC class II 


4 


2.20E-02 


Cellular 
Component 


MHC class II protein 
complex 


6 


3.70E-04 


Cellular 
Component 


MHC protein complex 


7 


5.00E-03 


Cellular 
Component 


large ribosomal subunit 


3 


2.90E-02 


Molecular 
Function 


transcription regulator 
activity 


20 


1 .80E-03 


Molecular 
Function 


sequence-specific 
DNA binding 


13 


2.60E-03 


Molecular 
Function 


ligand-dependent nuclear 
receptor activity 


7 


5.80E-03 


Molecular 
Function 


transcription factor activity 


15 


1 .20E-02 


Molecular 
Function 


steroid hormone 
receptor activity 


6 


2.10E-02 


Molecular 
Function 


DNA binding 


20 


2.40E-02 


Molecular 
Function 


phosphatase regulator 
activity 


4 


2.90E-02 


Molecular 
Function 


protein phosphatase 
regulator activity 


4 


2.90E-02 



these miRNAs mentioned above may participate in devel- 
opment of the piglet digestive tract. 

Notably, some miRNAs among the top 10 identified 
here have been reported to be related to immunity 
(miR-320, miR-181a, miR-30a-3p, let-7a, let-7f and 
let-7c) and development (miR-193a-3p, miR-378 and 
miR-191). MiR-193a-3p was demonstrated to regulate cell 
proliferation, cell cycle progression in vitro and in nude 
mice [69]. MiR-378 promotes osteoblast differentiation by 
targeting polypeptide N-acetylgalactosaminyltransferase 7 
{GalNAc-T7ov GalNTT) [70], and miR-191 regulates eryth- 
roid differentiation in mammals by up-regulating erythroid- 
enriched genes RiokS and Mxil [71]. Meanwhile, miR-320 
is able to inhibit HL-60 cell proliferation by suppressing 
receptor 1 {TfR-1; CD71) [72], and miR-lSla was believed 
to act as an intrinsic antigen sensitivity "rheostat" 



Table 8 KEGG pathway analysis of potential targets of 
top10 miRNAs 



Term Count gene p-value 



cci^nA'^IA- (~ oil ^Hhiocir^n 
bbLU'4-J 1^. ^^t:ll dUMcblUII 

molecules (CAMs) 


1 4 


FUR, LOCI 00521 555, 
SELE, SELL, SELP, SLA, 
SLA-DMA, SLA-DOA, 
SLA-DOB, SLA-DRA, 
SLA-DRBl 


1 .OJC Uj 


SSC05310: Asthma 


7 


CD40, SLA, SLA-DMA, 
SLA-DOA, SLA-DOB, 
SLA-DRA, SLA-DRBl 


5.55E-03 


SSC04672: Intestinal immune 
network for IgA production 


8 


CD40, SLA, SLA-DMA, 
SLA-DOA, SLA-DOB, 
SLA-DRA, SLA-DRBl, 
TGFB3 


2.15 E-02 


SSC04330: Notch signaling 
pathway 


5 


DLL4, PTCRA, 
LOC733643, 
APHIA, N0TCH4 


2.33E-02 


SSC04612: Antigen processing 
and presentation 


10 


CD4, LOCI 001 52370, 
NFYB, PSMEl, SLA, 
SLA-DMA, SLA-DOA, 
SLA-DOB, SLA-DRA, 
SLA-DRBl 


3.00E-02 


ssc00564:Glycerophospholipid 
metabolism 


6 


AGPATl, AGPAT4, 
IAGPAT6, GNPAT, 
LOCI 001 52491, 
PCYTIB 


4.25 E-02 


ssc00020:Citrate cycle 
(TCA cycle) 


5 


AC02,DLST, IDH2, 
LOCI 001 57889, 
PCKl 


4.57E-02 



during T cell development [73]. MiR-320, miR-181a, 
miR-30a-3p and let-7 were shown to be downregu- 
lated in colorectal cancer [74], Of course, further ex- 
perimental evidence is needed to verify that these 
miRNAs are indeed related to immunity of the piglet 
digestive tract. 

IgA is a major immunoglobulin in milk [75]. Expres- 
sion of the polymeric IgA receptor (pIgR) in mammary 
epithelial cells contribute much to the development of 
the immune system at the early stage of lactation [76]. 
In the present study, some miRNAs were predicted to 
target genes {CD40, SLA, SLA-DMA, SLA-DOA, SLA- 
DOB, SLA-DRA, SLA-DRBl and TGFB3) involved in 
processes of the intestinal immune network for IgA 
production in porcine milk. CD40 is a B-cell antigen 
activated during immune responses [77]. CD40 and 
CD40 ligand {CD40L) expressed on activated T cells 
are essential to B cell proliferation [78] and secretion 
of IgG, IgA and IgE [79]. SLA Class I were found to 
be expressed in the epithelial and lamina propria 
cells of the intestine in adult pigs and to be in- 
volved in mother-newborn interactions [80]. A study 
in humans showed that TGF-^ acts as a specific switch 
for IgA present at early stages of development of B 
cells [81]. 
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Figure 13 MiRNAs targeting tiie IgA immune networlc. The top 20 miRNAs were analyzed, and 14 of them were found to participate in the 
IgA immune networl<, involving 20 target genes. Different colors and shapes represent various relationships between miRNAs and genes. 



In the present study, the top 20 miRNAs were used 
for IgA network analysis. APRIL was the predicted 
target of miR-193a-5p, which is essential to triggering 
IgA2 class switch in human B cells. Intestinal epithe- 
lial cells (lECs) release APRIL after sensing bacteria 
through Toll-like receptors, and mucosal vaccines acti- 
vate lECs to induce more effective IgA2 responses 
[82]. The let? family and miR-423-5p were predicted 
to target CCL25, a potent and selective chemoattract- 
ant for IgA antibody-secreting cells [83]. CCL25 is 
known to selectively modulate immune responses, 
specifically the localization of T lymphocytes to 
the small-intestinal mucosa [84]. CD80 and CD86, 
which are costimulators of T lymphocytes [85], were 
identified as possible targets of five miRNAs in our 
study. Let-7a, let-7c, miR-181b, miR-185, miR-378 and 
miR-423-5p were predicted to target the inducible co- 
stimulatory molecule (ICOS), which plays a key role 
in regulating T-cell differentiation, T-cell proliferation, 
and secretion of lymphokines, providing effective help 
for antibody secretion by B cells [86]. We hypothesize 
that some miRNAs identified here in porcine milk 
regulate IgA production in the intestine of piglets, 
which may play an important role in mucosa immun- 
ity. However, their regulatory mechanisms warrant fur- 
ther study. 



Conclusions 

In conclusion, the present study revealed 176 known 
miRNAs and 366 (315 mature miRNAs) novel pre-miRNAs 
in porcine milk, most of which were predicted to be involved 
in regulation of digestive tract development and immunity of 
newborn piglets. These findings contribute to an increased 
understanding of the roles of miRNAs in porcine (S. scrofa) 
milk exosomes and to building the foundation for under- 
standing their physiological functions and regulatory 
mechanisms. 
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